THE EFFECTIVE BANDWIDTH PROBLEM REVISITED 
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f^ . Abstract. The paper studies a single-server queueing system with au- 

tonomous service and £ priority classes. Arrival and departure processes 
are governed by marked point processes. There are I buffers correspond- 
ing to priority classes, and upon arrival a unit of the fcth priority class 

^ I ' occupies a place in the fcth buffer. Let A'^''"', k = 1,2, ... ,1 denote the 

quota for the total fcth buffer content. The values A'*"' are assumed 
to be large, and queueing systems both with finite and infinite buffers 
are studied. In the case of a system with finite buffers, the values A'*^' 

^ ] , characterize buffer capacities. The paper discusses a circle of problems 

related to optimization of performance measures associated with over- 
flowing the quota of buffer contents in particular buffers models. Our 
C^ . approach to this problem is new, and the presentation of our results is 

simple and clear for real applications. 
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2 vyacheslav m. abramov 

1. Introduction 

1.1. Approach. During the last two decades there has been an increasing 
interest in the effective bandwidth problem for queueing systems with prior- 
ities. There are different classes of messages (units) arriving in telecommu- 
nication systems, and all of them are characterized by their quality of service 
requirements. In order to provide these quality of service guarantees and to 
allocate necessary network resources, different priority classes characterizing 
units arriving to that network are used. 

There are a large number of papers related to this subject. A detailed 
review of the related literature (up to publication time) can be found in 
Berger and Whitt [9] (for further discussions see also [TO] ) . For other relevant 
contributions to this subject see [TT], [TO], f22], [23], \M\, [35], [38]. 

These papers all discuss approximations and suggest algorithms for opti- 
mal solutions for the allocation of resources or effective bandwidth problems. 
Most of these papers use large deviation techniques. For example, Elwalid 
and Mitra [21], [22] use Chernoff's inequality to approximate loss probabil- 
ities in finite buffer systems with large buffers. Berger and Whitt [9] also 
use exponential asymptotics [TO], [l5] for the workload high level crossing 
of the ith class priority unit. Other papers (e.g. [11], [32], [37], [H], [39] ) 
also apply one or other techniques of the large deviation principle. Many of 
the aforementioned papers are aimed at solving concrete analytic problems, 
and their results are based on an analysis of analytic transformations (such 
as Laplace-Stieltjes or the z-transform) and their approximations. Many 
of these results are then applied to M/G/l-oriented queueing models or to 
models with more general arrival processes having a Markov structure. 

The approach of the present paper substantially differs from these previ- 
ous ones. The main focus of this paper is the solution of bandwidth problems 
for G//M/l-related priority systems. To the knowledge of the author, such 
priority systems are not presented in the literature where the overwhelming 
majority of priority queueing systems studied are of M/GI/1 type. The 
innovations of the present paper are as follows. 

1. We consider models of queues with an autonomous service mecha- 
nism (see e.g. [TO], [TO] as well as Section [L3l of this paper). The main 
results of our analysis are based on stochastic equations, and our models are 
studied under a rather general setting and can be applied to a broad class 
of real telecommunication systems. The obtained stochastic equations are 
then used for analysis of particular systems with exponentially distributed 
service times, which are a subclass of queues with an autonomous service 
mechanism. Note that martingale techniques for priority queueing systems 
(different from the systems considered here) have been developed by Kella 
[30] . However, the approach of [30] differs from the present one. Specifically, 
[30 ] studies fluid networks of parallel queues with dependent Levy inputs. It 
shows that the special construction given in the paper can be applied to the 
analysis of workload processes in M/G/1 queues with a preemptive resume 
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discipline. The paper of Kella [30] is based on an extension of the earher 
results of Kella and Whitt [31]. In contrast, our approach is based on a di- 
rect construction of queues with autonomous service mechanisms, and can 
be applied both to M/G/1 and GI jMjX oriented priority queueing systems. 
[GI jMjX queues are precisely described in the paper.) 

2. The buffer content process is described by the so-called buffer type 
stochastic equation. The buffer type stochastic differential equation is a 
special stochastic differential equation with discontinuous right-hand side 
(see Filippov [24J) and has already been used by Elwalid and Mitra ^21j . 
[22] to study the model with two priority classes. However |21j and [22] 
used the explicit forms of this equation related to low and high priority 
units. Analysis of these explicit equations is a hard problem. In contrast, 
our buffer type equation are represented in an (equivalent) integral form, 
and we discover a very simple representation for cumulative buffer contents, 
see Lemma 12.11 and Theorem 12.21 According to this representation, the 
system of equations for cumulative buffer content processes is the usual 
system of stochastic equations describing standard queue-length processes 
with an autonomous service mechanism. This finding essentially simplify the 
analysis, algorithms of solution and finally gives very simple approximation 
of the explicit solution. For example, it enables us to study the system with 
an arbitrary number of priorities. 

3. Some papers (e.g. [21], [22]) assume that buffers have large capacities 
and discuss the probabilities of buffer overflow. They use general estimates 
given by large deviation theory, and particularly, by Chernoff 's inequality. 
Being well-motivated theoretically, these estimates do not properly solve real 
practical problems. There is an example in [18] showing that inequalities 
based on exponential bounds can give unrealistic results. 

We offer a unified approach to systems with finite and infinite buffers. 
Large parameters N^'^', k = 1,2, ...,£, that are used in the sequel, are 
referred to as quota for buffer content and are related to finite and infinite 
buffers systems. In the case of finite buffers models with recurrent input 
and exponentially distributed service times of batches, we develop the known 
asymptotic results on losses in GI/M/l/n queues as n — > oo [2] to the case of 
GI/M^='^/l/Nk queues (A; = 1, 2, . . . , ^) with large buffers N^ (the second 
position M of the notation GI/M /l/N^ means that the service 

time of units is exponentially distributed, and batch size is equal to C) and 
then adapt the obtained asymptotic result to estimate the loss probability 
in systems with large finite buffers. The asymptotic representation of this 
paper, that is used for the loss probability in GI/M^^'-^ /l/N^ queues and 
then for the probability of buffer overfiow, is preferable to general type 
estimates such as Cramer or Chernoff inequalities. The asymptotic results 
for the loss probability in GI/M /l/N^ queues are expressed via the 
roots of the appropriate functional equations (see Sections [5] and [6]) . They 
are also useful in studying the behaviour of losses in the case of heavy 
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load conditions. Cramer and Chernoff inequalities are rougher, but their 
advantage is that they are explicit. However there are exact estimates in the 
form of explicit inequalities in the literature for the stationary probabilities 
of GI /M/l/N large buffer queueing systems as well (see, |I6|), and they 
can be easily adapted to the loss probabilities of the standard GI /M/l/n 
queueing system (see the discussion section in [2] ) with application to models 
such as the GI /M^^^ /1/Nj^ queues considered in the paper. We however 
are not going so far. 

The results of this paper can be also applied to M/G^^'-' /\/N oriented 
large buffers models. However, in this case a special asymptotic analysis 
similar to that given in Abramov |i] is necessary. This asymptotic analysis 
is routine and not provided in the paper. 



1.2. Convention on the notation. For any increasing random sequence 
of points ti, t2r • • , the associated point process Z{t) = Yl^=i'^{^n < t\ is 
always denoted by a capital Latin letter. If Cij C2, ... is a sequence of marks, 
then the associated marked point process Z{t) = Ylij=i Cj is always denoted 
by calligraphic letters. All processes considered in the paper are assumed to 
be right-continuous having left-limits and starting at zero. Exceptions from 
this rule are especially mentioned in the text (e.g. Remark II. ip . For an 
arbitrary point process Z{t), its jump in point t is denoted AZ{t) = Z{t) — 
Z{t—), where Z{t—) is the left-limit of the process in point t. For arrival 
processes we use letters A and A with sub- or super-script (the notation 
is given in Section II. 3p . and for departure process we use letters D and 
T>. The buffer processes describing the buffer contents will be denoted by 
calligraphic letter Q with sub- or super-script (the notation is in Section 
II. 3p . All processes of this paper are assumed to be given on a common 
filtered probability space {il,T,F = {J^t)t>o,^}- 



1.3. Description of the system. The paper is concerned with priority 
queueing system having I buffers. Units arrive at the kth buffer at random 
time instants t\ = t^^ , t2 = t^ + T2 , ■ ■ ■ , and the nth unit arriving 

(k) 

at the kth buffer has a positive integer random length ■!?„ . (In telecom- 
munication systems length can represent required memory for the message.) 

Denote A^^>{t) = X^^^ I -^ in ^ t>, where !{■} denotes an indicator of the 

event, and A^'^^t) = ^fL^^ ^f- 

The departure process I'(i) is assumed to be a point process with constant 
positive integer jumps C. Let xi) X2,- ■ .denote times between departures, 
and let x„ = Y17=i Xi denote the nth departure moment. Then 2?(i) = 
^ 'l2'^=i^{^n < t}- The constant C is called depletion rate. 

The buffers are numbered 1,2,. . . ,£, and the buffer with lower order num- 
ber has higher priority. Assume that the buffers are infinite. Then the 
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equation for the first buffer content {highest priority buffer) is 

(1.1) g(i)(t) = maxjo, Q!-^\t-) + l\A^^\t) - AP(t)} . 

According to (jl.lh . the buffer content Q^'^'{t) is governed by the processes 
A^^\t) and T>{t) and is referred to as a queueing process with an autonomous 
service mechanism. Queues with autonomous service mechanism were intro- 
duced and originally studied by Borovkov [12], [13]. For different applica- 
tions see P, [3], [5], [7], [25], [26j and [27]. The term /\A'^^\t) is called the 
arrival jump at time i, and the term AX'(t) is called the possible departure 
jump at time t. The prefix possible underlines the fact that departures can 
occur only if the system is not empty. For further simplifications, through- 
out the paper we assume that arrival and departure processes are disjoint., 
i.e. the probability of simultaneous arrival and departure is 0. 

If t is a jump point of the process I'(t), then the real departure jump at 
time t is min{Q''-'^)(t— ), C}. Thus, if Qy^'{t—)={), then there is no departure 
jump. 

Q}^>{t) is the second buffer content, the priority of which is lower than 
that of the first buffer. Q^'^' (t) satisfies the equation: 

(1.2) Q^^\t) = maxjo, Q^^\t-) + A^^^)^^) 

- AP(t) - Q«(t-)]l{QW(i) = 0}}. 



Despite the fact that equation (II. 2p has a more complicated form than equa- 
tion (jl.ip . both of these equations are of the same type. The term /\A^'^'{t) 
is an arrival jump at time t. The structure of the departure jump is more 
difficult. For simplicity we discuss the case (. = 2 below. If t is at a depar- 
ture jump and {Q^^'{t) > 0}, then {Q^^'{t—) > C}, and the jump is related 
to the first buffer only. Otherwise, if {Q'^^(t) = 0}, then the following two 
cases are possible: 

(i) {0<Q«(t-)<c}, 

(ii) {Qa)(t-) = o}. 

In case (i) departures occur from the first buffer, the first buffer is completely 
emptied, and if the second buffer is not empty, then in the case {Q^^\t—) < 
C} departures occur also from the second buffer. In case (ii) departures 
occur merely from the second buffer, provided that this buffer is not empty. 
Thus the real departure jump in this case is 



mm 



{Q«(t-) + Q(^)(t-),c}. 
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Equation ()1.2p is easily extended to the kth buffer content for any k 
=1,2,. . . ,i. Indeed, denoting 

(1.3) Qk{t) = Q«(t) + Q(2)(t) + ... + g(fc)(t), 

(1.4) Ak{t) = A^'\t)+A^^\t) + ... + A^^\t), 
we have the following equation (k = 1,2, . . . ,i — 1): 

(1.5) Q^'^+i) (t) = max {o, Q^'^+i) (t-) + A^^'^+i) [t] 

- AV{t) - Qk{t-)\l{Qk{t) = 0}]. 

The extension of (|1.2p given by (jl.Sp is quite clear. The term AA^^~^^'(t) 
is an arrival jump at time t (if any) to the buffer content Q^^^^'{t—). The 
other term of ( 11.51) 

^AV{t)-Qk{t-)\l{Qk{t) = 0} 



is also similar to the corresponding term of (11. 2p . If t is a jump point, then 
the meaning of Qkit—) is the total content of all buffers, the priority of 
which is greater than the priority of the given k + 1st buffer before the jump 
at point t, and {Qk{t) = 0} is the event, that all buffers, the priority of 
which is greater than the priority of the given A: + 1st buffer, are empty after 
the jump at time t. 

In the sequel the process Qk{t) is called the A:th cumulative buffer content. 

1.4. Formulation of the problems. The paper is concerned with the fol- 
lowing problems. Let N^^' , N^"^' , ...N^^> be large positive integer values. 
Assuming that appropriate limits in probability exist, denote 

A('=)(t) 

A(^{t) is the total number of arrivals until time t. Then J^^' is the fraction 
of arrival instants when the length N^ ' of the fcth buffer is exceeded. Let 
a^^\ a^^v . . , a^' be real positive numbers, denoting cost rates, and 
(1.7) J = a«J« + a(2)j(2) + ... + „(^)jW. 

Typical questions arising here are the following. 

1. Assume that the parameters N^^' , N^"^' , . . . , N^^' are given, but the 
depletion rate C can be controlled. Under what value of the depletion rate 
C we have J < e, where e is a given positive small value? This question can 
be formally written as follows: minimize C subject to J < e. 

2. Assume that C is given, but A^*-^', A^^-^', . . . , A"'^-* are control variables. 
Assume additionally that with given l3^'^\ P^^\- ■ ■ i/?*-^-* the values A^^-"^), A^(^), 
. . . , A^(^) must satisfy the condition: A^^^) = [/3(2)A^(2)j = L/3(3)Ar(3)j = . . . = 
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\_P^ >N^ ' \, where [-J is the notation for the integer part of number. The 
problem is to minimize A^^^^ subject to J < e. 

Remark 1.1. (jl.6p apphes to the finite and infinite buffers systems. To finite 
buffers model we prescribe that a complete arrival group is rejected when 
upon arrival the buffer overflows. In the case of the system with infinite 
buffers, Q^'''{t), A: = 1, 2, . . . , ^ all are assumed to be right continuous having 
left limits. In the case of finite buffers model, Q^ '{t), k = 1,2, ... ,£ are 
not longer right-continuous. For example, if N^^' is the capacity of the 
first buffer, and at moment t- the buffer overflows, then we admit that 

Q'^^ itj 1 is greater than A^'^^ in t- . However in the neighborhood of this 
point Q(^) (tf^^ < iV(^). Then the left and right hmits of Q^^^*) in Point 
t- are not greater than N^^^' , both these limits (with probability 1) are 
equal and t^- ' is an isolated point. 

1.5. Brief description of the mathematical ideas, methodology and 
contribution of the paper. In this section we describe the mathematical 
ideas of this paper, as well as the methodology and overall contribution. 

We start from the description of the buffer content process. For the 
highest priority buffer the equation for the buffer content is very simple. It 
is described by equation (jl.ip . The equations for the lower priority buffer 
contents are relatively more complicated and described by equation (jl.2p . 
However, the equations for the cumulated buffer contents are simple and 
described by a difference recurrence equation similar to (jl.ip 

(1.8) Qk{t) = max{0, Qk{t-) + AAk{t) - AV{t)}. 
Another form for (jl.Sp is a stochastic equation 

(1.9) Qkit) = Ak{t) - V / I{Qk{u-) > j}dD{u). 

(In all these two equations k = 1,2, ... ,i.) The stochastic equation (|1.9p 
can be rewritten 

Quit) = Akit) - V{t) + V / I{Qfe(n-) < 3 - l}dl)(u), 

with subsequent reduction to a Skorokhod problem (see [8], [40], [H]). (In 
the case C = 1 such reduction was provided in [33] . For its further applica- 
tion see also [l] and [3].) 

It is shown then that representations similar to (jl.Sp and (jl.Op remain 
valid for finite buffer models. Thus, in all cases the problem reduces to 
analyzing queueing systems with an autonomous service mechanism. 

We use these results for analysis of particular queueing buffer models 
with priorities. A system with exponentially distributed service times is a 
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special case of a system with an autonomous service mechanism. (A special 
construction of models with finite and infinite buffers is explained later in 
Section [5] of the paper.) 

In the case of finite buffer models with renewal input and exponentially 
distributed service times we adapt recent results on asymptotic analysis 
pj. As in papers [2], [4] and [6], the analysis is based on reducing the loss 
probability to a convolution type recurrence relation: 

n 
fn = 22 f-ri-j+l'^j (/o > 0), 

i=o 

where ttq > 0, ttj > for all j = 1,2,..., and X^^o ^i ~ ^' ^^^ applying 
asymptotic analysis similar to that of the book of Takacs [43j, p. 22-23. 
Consequently, we provide heavy traffic analysis of these models based on 
asymptotic expansions of the results obtained under "usual" conditions. The 
loss probability for the large finite fcth cumulative buffer is then not greater 
than the sum of the loss probabilities in the associated GI jM^^^ jXjNi 
queues, z = l,2,...,A;. However, for large values Ni this sum is very small, 
with the order of this sum being the same as the order of one (maximum) 
term obtained by asymptotic analysis, and an estimate obtained seems to 
be better than that estimate obtained by rough methods of large deviation 
principle and Chernoff 's inequality. 

Thus, the main mathematical contribution is a general theory of priority 
buffer models with application to particular priority queueing systems with 
recurrent input and large buffers. 

1.6. Organization of the paper. In Section [21 Lemma [2T] states that the 
fcth cumulative buffer content has the representation (j2.ip . The intuitive 
sense of Lemma [2. II is that the fcth cumulative buffer content for the system 
with infinite buffers is described by the same equation as the queue-length 
process in the queueing system with autonomous service mechanism, the 
arrival process of which is A^^it) and the departure process I^(t). We fur- 
ther prove a stability theorem. The main condition for stability is (j2.4p . the 
proof being based on reduction to the Skorokhod refiection principle and 
results of Borovkov [12], [13]. In section [3] the model with finite buffers is 
considered. It is shown that the equations for cumulative buffer contents in 
this case are similar to the case of a model with infinite buffers. In Section 
m we derive the formula for J^^'^ k = 1,2,...,£, using the level-crossing 
method based on representation (j4.3p . In Sections [5] and [6] special models 
of queueing systems are studied. The results of these sections are illustra- 
tive, and we do not discuss general buffer models with batch arrival such 
as GI'^ /M^^^ /I queues, although the asymptotic geometrical bounds for 
stationary probability to reach high level N in GI /M /I queues is known 
(see [20]). All models considered here are particular cases of the general 
models discussed in Sections [2] and O these models are with independent 
identically distributed inter arrival times. The results of Section [6] are based 
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on an extension of recent results [2j. As in ^ the asymptotic analysis is 
based on reduction to appropriate representation helping us to use then the 
Takacs theorem on asymptotic behavior of the convolution type recurrence 
relation [43], p. 22-23. In Section [6.11 the asymptotic behaviour of losses are 
studied under "usual" conditions, while in Section 16.21 the analysis of losses 
is done under heavy load conditions. In Section [7] approximation of the ini- 
tial problem stated in Section [1.41 by another related problem is suggested. 
In Section [8] algorithms for numerical solution of the problems of Section [7| 
are proposed. There are concluding remarks in Section [9l 

2. The stability theorem for the infinite buffers system 

The representation for the buffer content of infinite buffers systems given 
by (jl.ip . (|1.2p and (|1.5p is difficult to analyze. However, for the cumulative 
buffer contents of infinite buffers systems the representation is simple. 

Lemma 2.1. For all k = 1,2, ... ,i the following equation for the kth cu- 
mulative buffer content Qk{t) holds: 

(2.1) Qk{t) = Akit) - V / I{Qk{u-) > jjdDiu), 

,=1^0 

where D{t) = -^■ 

The proof of this lemma is given in Appendix A. 

The statement of Lemma 12.11 has a simple intuitive explanation. For 
example, in the case i = 2 we have two classes of units, and clearly the 
cumulative buffer content process Q2(t) = Q^^'{t) + Q^'^'{t) contains two 
unit classes together, and therefore must behave as a usual (i.e. without 
priorities) queue-length process with an autonomous service mechanism, the 
arrival process of which is .42 (t) = A^^'{t) + A^'^'{t), and the departure 
process is 'D(t). This intuitive explanation is easily extended to the case of 
arbitrary k = 1,2, ... ,i number of classes. 

The right-hand side of this equation contains the sum 

V / I{Qk{u-)>j}dD{u). 

Nevertheless, the problem can be reduced to the Skorokhod reflection prin- 
ciple. 

Denote Sk{t) = Akit) - V{t), k = l,2,...,i. Then, 

(2.2) Qkit) = Skit) + V / I{Qkiu-) < j - l}dDiu). 

Equation (12. 2p implies that Qkit) is the normal reflection of the process Skit) 
(5fc(0) = 0) at zero. More accurately, Qkit) is the nonnegative solution of 
the Skorokhod problem of the normal reflection of the process Skit) at zero 
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(see Skorokhod [lO] as well as Tanaka [Sj and Anulova and Liptser ^, 
Ramanan 1381). This is because the function 



<t'k{t) = i2 [ HQkiu-) < j - l}dD{u) 

satisfies the following two properties: 

(a) Jq h[Qk{u)]d(j)k{t) = for any continuous nonnegative function h{x) 
with /i(0) = 0; 

(b) the function Jq[Y{u) — Qk{u)]d(j)k{u) is not decreasing for any non- 
negative right-continuous function Y{u) having the left limits. 

Let us show (a). We have 

t f-t 

h[Qkiu)]dcl)k{u) = / h[AQkiu)]d(l)kiu) 
Jo 



l\[AQk{u)]d Iy. f ^i^kiv-) < j - l}dD{t 



Let Ui denote the points of jump of the process Qk{u) in the interval [0,t]. 
For the last integral we have the following representation: 

I' h[AQk{u)]d I J2 fHQkiv-) < j - l}dD{i 
c 

0<Ui<t j = l 

The last sum is a finite sum: the number of points Ui is finite in any finite 
interval [0,t] with probability 1. Any value of jump AQfc(nj) is bounded 
with probability 1, and the nonnegative continuous function /i[AQfc(u)], 
satisfying the property /i(0) = is therefore bounded for all < n < t. 
In addition, taking into account that the jumps of the processes Ak{u) and 
D{u) are disjoint, i.e. either AAk{ui) = or AD{ui) = with probability 
1, we arrive at the conclusion that J^ h[Qk{u)]d<j)k{u) = 0. (a) follows. 

(b) is implied by (a). 

It follows from the Skorokhod reflection principle that the function 0fc(t) 
has the following representation: 

(pk{t) = -iniSkiu). 

u<t 

Therefore Qk{t) has the following representation 

Qk{t) = Sk{t)-iniSk{u), 
(2.3) "^* 

A: = l,2 t 
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Equation ()2.3p is well-known in queueing theory. Following Borovkov [12j , 
we have the following statement of the stability. 



Theorem 2.2. Assume 
(2.4) ] 



,. Ai(t)-V(t) 

hm —^--^ ^— = r < 

^t^oo t 

and Sk{t), k = l,2,...,i, are stationary point processes, the increments 
of which coincide in distribution with the corresponding increments of the 
processes Ak{t) — V{t), k = 1,2, . . . ,i. 

Then there exist stationary processes Q^^'{T), k = 1,2,. . . ,£, such that 



(2.5) 
and 

(2.6) 



Q«(r)^sup Si{T)-Si{u) 

u<T L 



Q(^)(r)=sup Sk{T)-Sk{u) 

u<T L 

rC Z^ O, . . . 5 -C. 



sup 

u<T 



Sk-i{T) - Sk-i{u) 



Proof. The proof is based on representation (|2.3|) and can be found in 
Borovkov [12]. Specifically, it follows from that proof that there are sta- 
tionary processes Qk(T), k = 1,2, . . . ,i such that 



(2.7) 



Qk{T)= sup Sk(T) - Sk{u) 



u<T 



Therefore, keeping in mind that Q^'^^t) = Qk{t) - Qk-i{t), k = 2,^,...,i 
and Q^^'^it) = Qi{t), from 1^ we have ([22]) and i^Ml- □ 



3. The finite buffers model 

Equation (j2.ip and other related equations for infinite buffers content can 
be easily extended for the model with finite buffers. It is assumed that 
if upon arrival of a batch the buffer of a given class overflows, then the 
complete arrival batch is rejected, see Remark 11.11 

For the analysis of the finite buffers case we introduce new arrival pro- 

cesses A (t), which are derived from the initial processes A^ '{t) as follows. 
We set 



(3.1) 



AA^''\t) = A^(*^)(t)I |QW(t) < N^^^l . 



-ii(k), 



The arrival processes A (t) take into account only jumps of real buffer 
content process. Thus A{t) — A{t) is the number of lost units during time t, 
and A{t) — A{t) is their total length during that time t. 

Then the buffer content process Q^^^ (i) is defined by the pair of equations 

(3.2) Q^'\t) = maxjo, Q^'\t-) + AA^'\t) - AP(t)} , 
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(3.3) Q^^\t+) = maxjo, Q^^Ht-) + AA^^\t) - AP(t)} . 

Thus in the case A^(i)(t) = AA^^\t) the buffer contents Q'^^\t) and 
Q^^'{t+) are equal and there is no loss at time t. Otherwise, if AA^^'{t) ^ 

AA^^\t), i.e. AA^^\t) = and AA^^\t) > 0, then there is a loss of a unit 
in time t. 

Next, similarly to (jl.5p for k = 1,2, . . . ,i — 1 we have another pair of 
equations: 

(3.4) Q^^^+i) (t) = max {o, Q^^'+i^ (t-) + A^^'^+i) (t) 

- AP(t) - Qk{t-) - AAk{t)]l{Qk{t) = 0}}, 

(3.5) Q('=+i)(t+) = max{o,Q('=+i)(t-) + A:4^''+'^(i) 

- AP(t) - Qk{t-) - AAk{t)]l{Qk{t) = 0}}. 

Similarly to Lemma [2. H for the finite buffers model we have the following 
lemma. 

Lemma 3.1. For all continuity points of the kth cumulative buffer content 
process Qk{t), k = 1,2, . . . ,£, we have: 

(3.6) Qk{t) = Mt) - E / HQkiu-) > j}dD{u). 

4. The formula for J*^'^^ 

In this section we study the dynamics of the buffer lengths by level- 
crossings analysis for the infinite buffers model. It is assumed throughout 
that condition (j2.4p for the stability is fulfilled. 

In addition to the stability condition assume: 

(4.1) p/iii^^i^ = A('=)| = l, 

and 

(4.2) p/lim :^=/il = 1. 

t^co t 



1 ^^^?^ „q(A:) 



Then according to (j2.4p the sequences -J27=i''^i ' ^ ~ 1)2,..., £, 



as 



n — > oo, also converges with probability 1. 

Recall that t^ = rf \ 4*"^ = 4''^ + 4''\ . . . , {k = 1,2, . . . ,i) denote the 
sequence of points (arrival moments) of the process A''''\t), and xi = xi; 
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X2 = Xi + X2, ■ ■ .denote the sequence of points (the moments of possible 
departure jumps) of D(t). 

Then, for the number of up- and down-crossings for m > 1 we have the 
following equation: 



j=i 
Dit) 



(4.3) 



^ I |m < Q^'^) (xj-) <m-l + c\ 

i=i 

C D{t) 

j:j:i{Q^'^{x,-)=m-i+i} 

1=1 3=1 



+ l{Q^^Ht)>mY 



where Q}-^'{&) = 0. Equation (|4.3p can be explained as follows. The left- 
hand side of the equation is the number of arrivals until time t, seeing 
before arrival the buffer content less than m and at the moment of arrival 
not smaller than m. This constitutes the number of up-crossings of the level 
m until time i, i.e. the number of instants where arrivals jump over the 
level m — \. The first term of the right-hand side describes the number of 
departure moments when immediately before departure the buffer content 
is between m and m + C — 1. (Then after the departure the buffer content 
is between max{0, ttz — C} and ttt. — 1, and this constitutes the number of 
down-crossings of the level m). The difference between the number of up- 
crossings and down-crossings of level m can be either 1 or 0, and the second 
term of the right-hand side compensates for this difference. 

Dividing the both sides of M.Sh by t, and letting t increase unboundedly, 
we obtain: 



A(*=)(t) 
^Inn Ie Y. l{Q^'^(f^) >"i, Q('^)(tf)-j <m; 

(4.4) *"" * '■=' 

^ ' . C D{t) 

= lim -E^^l|Q(^')(a;j-) = m-l + /|, 
1=1 i=i 
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and after elementary transformations (see Appendix B) we arrive at 

X(k) 1 



jik) 



A(i) + A(2) + . . . + A(^) A(fc) 



(4.5) 



lim -E /" V I ( QC^) (u-) > iV^'^) + /l dD{u) 
^°^t Jo fr^ ^ J 

1 1 ft 

— lim-EV / i|qW(u-) >7VW + /|dL'(u) 



A; — 1, 2, . . . , t, 
where Xe = A^^) + A^^) + . . . + A^. 

5. The buffers content distribution of GI /M^^^ /I queues 

5.1. Main result. We start this section with a representation for the buffer 
content processes in the case where the arrival processes A^ '{t), k = 1,2,. . . , 
a. all satisfy (j4.ip . and the process D{t) is Poisson. Assume also that ■(?„, = 1 
for all n. 
We have: 

iim -E / 



(5.1) 



Iim ^E / l{Qk{u-) =m- l}d^fc(n) 

:^oo t Jo 

I '^ ft 

/i Iim - V / "^{Qkiu) = m - 1 + l}du, 



k = 1,2 



"-^ -^^ 



However ()5.ip does not permit us to obtain explicit results for the stationary 
probabilities even in the case where the processes A^ '{t), k = 1,2, . . . ,i all 
are renewal processes. Moreover, in the case where all processes A^ '{t), 
k = 1,2, . . . ,i are renewal, Lemma [2. II is no longer useful in general, because 
stationary interarrival times to cumulative buffers are dependent in general, 
and the corresponding stationary arrival processes are not longer renewal. 

Therefore we consider the following special case of the general buffers 
model. Let A(t) be a point process of arrivals satisfying the condition 
Pjlimt^oo^ = A} = 1. Let vr^^), 7r(2),..., vr^ be positive probabili- 
ties, X^fc=i tt'" = 1, where vr'^ is a probability that an arriving customer 
belongs to the class k. Then the points processes A^''\t), k = 1,2, . . . ,£, 
are all thinnings of the original process A{t), and in the case where A(t) is 
a renewal process all the processes A^ \t), k = 1,2, ... ,£ are renewal pro- 
cesses as well with intensities A*^^ = Avr'^ correspondingly. Consequently, 
the processes Ak{t), k = 1,2, ... ,i, are renewal processes with intensities A^ 
= A(^) + A(2) +...+ A^"), and one can apply the theory to the GI/M^=^ /I 
and GI /M^^^ /1/Nk queues with large buffers N^. 
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By GI /M /I queue we mean a single-server queueing system with 
recurrent input and exponentially distributed service time of the constant 
size batch C. In the sequel we use the notation G//M*^/l for these queue- 
ing systems. GI /M^ /I queueing systems are particular systems with an 
autonomous service mechanism, and they are therefore described by buffer 
type stochastic differential equations or by one of the above equivalent forms 
of these equations. For these queueing systems therefore Lemma [2.11 remains 
true. Specifically, from this lemma one can conclude that the cumulative 
buffer content processes are described by the steady-state distributions of 
the usual queue-length processes of the GI /M*-^ /I queues. The stability 
condition for these queues is pn = —^ < 1. 

Using a standard method, the limiting and stationary probabilities for the 
cumulative buffer contents of GI /M^ /I queues are calculated as follows. Let 
ifcj denote the jth arrival moment to one of the first k buffers. Then for 
limiting and stationary probability we have the following. 

Theorem 5.1. For cumulative buffer contents Qk{t), k = 1,2, . . . ,i. 

(5.2) Pk^rn = lim HQkitkj-) = m} = <?r(l - ^fc), 

where (,k is the (unique) root of the functional equation 

(5.3) z = Bkip-fiz^) 

in the interval (0,1), and Bk{s) = L e~^^dBk{x) is the Laplace- Stieltjes 
transform of the stationary distribution of interarrival time Bk{x) to the 
first k buffers. 

Proof. The stationary probabilities of GI-^ /M^ /I queues can be found in 
Economou and Fakinos [20], and the statement for G//M'^/l queues can 
be deduced from their result |J However, the direct proof of the result for 
the GI /M^ /I queue is much simpler than that reduction from the afore- 
mentioned general result. Therefore below the direct proof of this theorem 
is provided. 

First of all notice, that according to (15. ip the state probabilities immedi- 
ately before arrival, Pk,m-, are 

1 /■* 

Pk,m = hm --E / l{Qk{u-) = m}dAk{u), 

t->oo Afct Jo 

and Pk,m = {I — z)z"^ for some z < 1. Let fm denote the number of up- 
(down-) crossing of level m during a busy period of GI/M /I queue (the 
number of cases where immediately before arrival there are m customers in 
the system). Then, by renewal arguments E/^ = z"^, and according to the 



The following additional condition is missed in the main statement of [20]: the common 
divisor of possible values X and Y must be equal to 1. 
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total expectation formula for any ?7i > 1 we have the following equation: 
(5.4) ^"^ = 1;/ e-'^^-fc^ z™-Mi?,(x), 

where Bk{x) is the probability distribution function of interarrival time. 
Therefore, from (15.41) we obtain 






z^) 



dBkix) 



and the statement of Theorem 15.11 follows. By standard method (see e.g. 
|42j . |28j ) one can prove that under the assumption p£ < 1 there exists a 
unique root of equation z = B^^fj, — jiz^) in the interval (0,1). D 

5.2. Particular case. We consider an M/M /I queueing system with in- 
finite buffers. This particular case is easily deduced from the statement 
of Theorem 15.11 Specifically, in the case of Poisson arrivals from (j5.3p we 
obtain the equation: 

Then, the constant ft must be the solution of equation 

(5.5) - = E-^ 

belonging to the interval (0,1). A similar result can be also found in [T5] for 
nodes of a network, the customers of which are served by random batches. 

6. Loss PROBABILITIES FOR CUMULATIVE BUFFERS 

In this section we discuss loss probabilities assuming that the kth cumula- 
tive buffer content has large capacity Ni^.. We study buffer loss probabilities 
under "usual" and heavy load conditions. By "usual" conditions we mean 
the case when the load parameter of the queueing system is fixed, while in 
the case of heavy load conditions the sequence of load parameters, associated 
with series of queueing systems, approaches 1. 

6.1. Loss probabilities under "usual" conditions. We use the notation 
GI /Ai^' /1/N]^ for the queueing systems with finite capacity N}^, similar to 
the notation used for the queueing systems with infinite capacity in the 
previous section. According to Lemma [3. II the cumulative buffer contents in 
continuity points of the process Qk{t) behave as usual GI /M^ /1/Nk queues. 
However, the behavior of the number of losses, the main characteristic of 
interest, is essentially different, that is the losses in GI /M /1/Nk queues 
are not equal to the losses in the corresponding cumulative buffers Qk(t). 
Specifically, the losses in GI /M^ /1/Nk queues occur only in the case in 
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which the buffer overflowed when the arriving customer met all waiting 
places busy. The losses in the cumulative buffers Qkit) can occur in many 
cases when one of specific buffers, say jth buffer, 1 < j < fc, has overflowed. 

However, in some cases when the values A^i < A'^2 < . . . < A^^ all are large, 
a correspondence between GI/M^ /l/N^ queues and finite buffers models, 
may give useful asymptotic results. 

Specifically, the loss probability of a customer arriving at one of the first 
k buffers is not greater than pi + P2 + ■ ■ ■ + Pk, where pi denotes the loss 
probability in the corresponding GI/M'-'' /l/Ni queueing system, the prob- 
ability distribution of interarrival time of which is Bi{x). All probabilities 
Pk: ^ ^ k < i are very small as is Nf^ large. They decrease geometrically fast 
(see Theorem 16.11 below) , and the finite sum of these probabilities seems to 
remain a good upper bound for the buffers loss probability. 

Theorem 6.1. The buffer contents loss probability is not greater than pi + 
P2 + ■ ■ ■ +Pk, k = 1,2,. . . ,£, where 

(l-p,)[l + C/.B^.(M-A^?F)kf 



Pk 



(6.1) (1 - pfc)(l + ?fc + . . . + ?fc^-^) -Pk[l + CfxB'j^ifi - fi<;^)]4 

^k I ' 



+ o(,f^^^ 



Afc 
Pk 



C/i' 
and ?fc is the (least) root of the functional equation 

z = Bk{p- pz^) 

in the interval (0,1). 

Proof. We consider the GI/M /1/A^fc queueing system. Following Miyazawa 
[55] . the loss probability for the GI /M^ /1/Nj^ queueing system is deter- 
mined by the formula 

1 
Pk 



where the generating function of vr^ j, j = 1,2, . . . is 

(6.2) Uk{z) = ^ Tik,jz^ = g- — , 

j=o Bkiji- pY{z))-z 

and Y{z) is the generating function of complete service batch. In the case 
of the GI /M*-^ /l/N^ queueing system Y{z) = z'~^ , and (|6.2p can be then 
rewritten as 

(6-3) nfc(z) = ^ — . 

Bkifi- fiz^) - z 

In the particular case of C = 1 the asymptotic behaviour of the loss proba- 
bility has been studied in [2] and [17]. In the case of p^ < 1 it was based on 
an application of the Takacs theorem [l3], p. 22-23. 
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In the case of C > 1 the scheme of the proof is similar. Expanding {1-z ) 
in the numerator of (16. 3p as 1 — z = (1 — z)(l + z + . . . + z ) we have 



(6.4) Uk{z) 



{1 - z){l + z + z^ + . . . + z'--^)BkifJ, - fJ.z 



^C\ 



BkifJ, - fiz'^) - z 
Therefore, the other generating function Iik{z) = j^Hkiz) is 



(6.5) 



^k{z) = ^7rk,iz' = ^ X^^fc.i I z' 

i=0 j=0 \j=0 

(1 + z + z2 + . . . + z'^-^)Bk{fi - nz^') 



BkiiJ,- HZ 



c\ 



and the loss probability is 
(6.6) 



Pk 



'"'k,Nk 



Our goal is therefore to find the asymptotic behaviour of T^k,Nk ^^ ^k -^ oo. 
The equation z = B}^{jjl — ^iz^) has exactly one solution "j^ in the interval 
(0,1). Furthermore, -Bfc(/U — ^z*-^) is the probability generating function of 
some integer random variable, i.e. 



and 



where 



Bk{ii-Hz'^) = R{z) = Y, 

i=o 

_ c-i 

Iik{z) = F{z)Y,z\ 

3 

R{z) 



rjz^, 



i=0 



F{z) = Y, Uz^ 



n=0 



R{z) - Z 



Therefore (see Takacs [l3]), the sequence {/n} satisfies the recurrence rela- 
tion 

n 
fn = Y fn-j+irj, h > 0. 
3=0 

Since <^k < 1) we correspondingly obtain 7 = J2'^=i''^^n > 1- According 
to formula (35) of [43j, p. 23 



(6.7) 



lim 

n— >oo 



/n 



[1 - F'i5)] 



/o 



1-7' 



where 6 is the least root of equation z = R{z) in the interval (0,1). 
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In our case pk = ^, and i^k = <J, and we have: 



l + ?fc + ... + ?^"^ 1 



Pfc - 1 



TTfc.Affe 



and the statement of the proposition follows from (16. 6p after some algebraic 
transformations. D 

6.2. Loss probabilities under heavy load conditions. The loss proba- 
bilities under heavy load conditions for GI/M/l/n queues have been recently 
studied in [2] and [36] . For the further development of these results see also 
[6] and [47] . In this specific case, the behaviour of the system under heavy 
load condition differs from the classic cases considered in these papers. 

We consider the case of heavy load conditions and assume that the load 
parameter p£ = -^ is close to 1. More specifically, we assume that pe = 
pe{S) <1 {6 is a, small parameter), and pe{6) approaches 1 from the left 
as 6 vanishes. Denote pej = p^ j'^ x^dBi{x) {p£^i = jjz-)- We have the 
following result. 

Theorem 6.2. Assume that pe{S) < 1, peiS) approaches 1 from the left, 
and 5Ni{5) ^ A > as (5 vanishes. Assume also that pe^2 = hniiJ^o P£,2(<^); 
and pifi{6) remains hounded as 5 vanishes. Then 

(6.8) « = VIll^.U +„(„]. 

Proof. Let us first derive the expansion for the least root of equation z = 
B^iyp — pz'^) under the assumption of the theorem. Clearly, the root of this 
equation approaches 1 as 6 vanishes. Therefore, using the Taylor expansion 
of Bi{p — pz^) as 5 vanishes, we obtain the following equation for ^^ 

, = 1 _ (1 _ 5)(1 _ ^) + (2)^^.2(1 - zf + 0(1 - zf. 
Ignoring the last term o(l — z)^ we have the quadratic equation, the solutions 



of which are z = 1 and z = 1 — 5/ 



(2)^2 



Therefore we obtain 



5 



(6.9) ^, = i_^__ + o(5). 

(2)^2 



Notice, that representation similar to (j6.9p for the root of equation z = 
Bi{p — pz) (particular case where C=l) has been obtained in Subhankulov 
m, P- 326. 
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Next, the asymptotic representation ior pk, k = 1,2, . . . ,i is given by 
6.ip . For k = i the main term of asymptotic expansion of 



is given by 

A \ 



(5exp 



(2)/'^,2, 



(see [2j, [6j for details of the proof), and according to (j6.9p the main term 
of the asymptotic expansion of 

1 + <;-^ + . . . + ?f -^ 

is given by 

c-±. 

Next notice, that the expansion for the term 

is 

1 — exp 



A 



(?)P^,2^ 



[l + o(l)]. 



Therefore, asymptotic relation (|6.8p follows. D 



7. Approximation of the solution in particular cases 

In this section we discuss the approximation of the solution for the prob- 
lem stated in Section [1.41 that is to minimize functional (jl.7p containing the 
terms a^^' J^^> associated with buffer contents Q^^'{t), k = 1,2, ... ,L 

However, in all particular cases above the explicit solutions were obtained 
for the cumulative buffer contents Qk{t), k = 1,2, . . . ,i, and the solution of 
the problem in the initial terms seems to be hard. Therefore we formulate 
and solve the problem in new terms. This solution of the new problem is 
then used to approximate the desired solution of the initial problem. 

Let us first introduce new functional Jk instead of the J^''' [k = 1,2,..., 
i), which were introduced in Section [1.41 

Namely, let A''i, A'^2;- • • ,]^£ denote large integer numbers, A''i < A'2 < . . . < 
Ne. We set 



Jk = P- lim -— V I {Qfc(tfcj) > Nk} , 



A: = 1,2, 
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Note first (see (I4.5|) ) that for J^, k = 1,2, ... ,i we have the following 
representation: 

Jfc = f lim^EV / I{Qk{u-)>Nk + l}dD{u), 

(7.1) Xft^oot fr^Jo 

K = 1, z, . . . , £, 

where A^ = A'^-'^) + A'^^) + . . . + X^''\ k = 1,2, . . . ,i. The proof of representation 
(j7.ip is similar to the proof of (j4.5p with minor difference in the notation. 
Replacing functional (jl.Tp by 

(7.2) J = aiJi + 02^2 + ••• + «£-'£, 

we have then the following problems similar to the problems formulated 
above in Section [L^ 

1. Assuming that A^i, N2,. ■ . ,Ni are known, minimize C subject to J < e. 

2. Assume that C is known, but A'^i, A'^2v • • ,Ne are unknown. Assume 
additionally that with given /?i, /32,. . . ,/?£-i the values A''i, N2,. ■ ■ , Ni must 
satisfy the condition: A^i = ll3iN2\ = [/32N3\=. . . =l(3e-iNe\. The problem 
is to minimize A^i subject to J < e. 

The values ai, 02, ■ ■ ■ , otf, N\, N2,. ■ ■ , Nn are unknown, and by approx- 
imation of the solution of the problem we hope to find a correspondence 
between the vectors (ai, a2, ■ ■ ■ , a^) and (a'^^, a^^\ . . . , a^') and between 
the vectors (iVi, N2,. . . , Ni) and (iV^^), N^'^\. . . , iV^^^) such that the so- 
lution of the initial problems formulated in Section 11.41 and the problems 
formulated in this section would be approximately the same. 

Consider first the queueing systems with infinite number of waiting places, 
say GI jM^ jl queues. 

Notice, that a\ = a^^' and A^i = N^^' . According to Theorem 15.11 the 
expected queue-length of the /cth cumulative buffer content immediately 
before arrival of a unit is 



(7.3) lim - V m / I{Qk{u) = m}du 

*^°° * ~i Jo 



m=i -^0 ^-* 



From (j7.3p we have the following. Put 

?i 



„ _ 1-ft _ '^1(1 -'^2) 

P2,l 



(7.4) "^'^ ^2 ?2(l-?i)^ 

l-?2 
P2,2 = 1-P2,1, 

and then 

(7.5) a2 = aip2,i + 0(^^2,2- 
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Similarly to (TTI]) and (173]) for fe = 1, 2, ...,£- 1 we set 

_ ^A:(l — ft+l) 
(7.6) <?A;+l(l-?fc) 

PA:+1,2 = 1 — Pk+1,1, 

and 

(7.7) Qfc+i = akPk+1,1 + a^''"^"^Vfc+i,2- 

Let us now express the correspondence between the vectors {N^^' , N^'^\. . . , 
TVW) and {Ni, N2,. . . , Ne). Let /?(2), /3(3),. . . , pW be such the real numbers 
that nW = [/3(2)Ar(2)j = L/?(3)iV(3)j = ...= L/3(')iVWj. 

Then for the purpose of approximation the values /32, Ps,- ■ ■ ,(3i are taken 
as 

(7.8) /3, = _^_^, A: = 2,3...,£, 

and iVi = L/?2iV2j = L/33A^3j = • • • = [(3iNi\. 

For the queueing model with large finite buffers, say GI/M /l/N^, k = 
1,2, ...,£, the approximation is similar. Specifically, since the buffers are 
large, approximation for (ai, 02,- ■ ■ , ai) can be given by (|7.3p - ()7.6p . The 
values /32, P3,. . . ,Pi are assumed to be taken by the same relation (I7.8[) . 

8. Minimization algorithms for the functional J 

In this section we discuss the problem of minimization of the functional 
J defined by (17.2p . Ji, J2,. ■ ■ , Je depends on parameters C, A^i, N2,. ■ ■ , A'^. 



8.1. C is known while A'^i is unkno'wn. Assume first that C is known, 
A'^i is unknown, and the problem is to find the value A'^i minimizing the 
functional J in the buffer models, where explicit representation for the state 
probabilities as well as for J^, k = 1,2, . . . ,1 are known. These models are 
considered in Sections [5] and O 

To be specific we refer to the models of infinite buffers of GI /M'-' /I 
queues. The algorithm has the following steps. 

• Step 1. Calculate <;fc, k = 1,2,. . . ,L Recall that <;k is the root of the 
functional equation (j5.3p in the interval (0,1). For each k it can be calculated 
by the fixed point method or by one of other well-known methods, say direct 
search method or gold section method (e.g. see [29j). 

• Step 2. We have I geometric distributions obtained in Step 1, and there- 
fore one can compute the corresponding values A'^i, N2,- ■ ■ ,N£ at which each 
of the tails of the geometric distributions multiplied to the corresponding 
coefficient a^, k = 1,2, ... ,i will be less than e (i.e. a^Jk < £)• 

• Step 3. By using the known coefficients (32, Ps, ■ . . , /?£ one can find 
the value iv{°"'^'". jv{°"'^'' is the maximum amongst all minimal values of 
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A''! under which akJk < e for ah fc = 1, 2, . . . , L Specificahy, we have the 
system: 

minA^i : aiJi{Ni) < e, 

minA^i : a2</2(-^i) < e, 



minA^i : a£Ji{Ni) < e, 

and jv{°""^'' is the maximum amongst £ obtained values of A'^i . 

• Step 4- By using the same known coefficients (32, Ps, • • • , /?£ one can 
find the value jV"^^^*^. N^pp^"^ is the minimum amongst all maximal values 
of A'"! under which a^Jk < | for all A; = 1, 2, ... , £. Specifically, we have the 
system: 

miniVi : aiJi(iVi) < |, 
minA''i : a2<^2(-^i) < Tj 



minA^i : aiJi{Ni) < -, 



and _/V"PP^'' is the maximum amongst i obtained values of A'^i. 

• Step 5. We solve the following integer programming problem: 

minimize iVi : Ni"""'^' < Ni < N:^^^''' 

subject to ai Ji + 02 ^2 + • • • + aiJi < e. 

8.2. A'^i, N2,. . . ,N£ are kno^vn w^hile C is unknow^n. In the case where 
Ni, N2,- ■ ■ , N£ all are known but C is unknown the algorithm of the problem 
solution is the following. 

• Step 1. From the stability condition find the lower (integer) bound for 
C: 

-ilower ■ I ^-( ^^ 



flower = min <^ C : -^ < 1 

• Step 2. Find^fc, A; = 1,2, . . . ,1 

• Step 3. Compute the functional J. 

• If J > e, then find a new value C and repeat steps 1-3. These procedure 
should be repeated more and more while J > e. Since the upper bound of 
C is unknown, the value C should be found according to the special search 
procedure offered by Rubalskii [39] . 

Rubalskii [39] proposed the minimization algorithm for a unimodal func- 
tion on an unbounded set. The optimal algorithm is an extension of the 
standard Fibonacci procedure. 
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9. Concluding remarks 

In this paper we studied queueing systems with priority classes and infi- 
nite and finite buffers. We derived general type equations for buffer content 
processes assuming that service mechanism is autonomous. The results of 
general theory were then applied to special queueing models with exponen- 
tially distributed service times. These queueing systems are a particular case 
of systems with an autonomous service mechanism. For the model having 
large buffers we derived an asymptotic result for the loss probability. We 
developed an algorithm for a solution of the problem numerically. 
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APPENDIX A: Proof of Lemma [O 

We start from equation (jl.ip . In order to write this equation in the cus- 
tomary form of a stochastic equation, we use the process D{t) = Yl^=i I{^n < 
t}. The jumps of the process D{t) are equal to 1, and according to the def- 
inition, we have CD{t) = T>{t) for all f > 0. Then due to the assumption 
that arrival and departure jumps are disjoint, (jl.ip can be rewritten 



c ,.f 
(Al) Q^^){t) = A^^\t)-Y,j i{Q^^\u-) > jYd 



,=i-o 



The equivalence of representations of (jl.ip and (A.l) can be easily checked 
by considering a small time interval {t — d, t] containing exactly one event as 
either arrival or departure of a unit. Then the term 

(A2) i{q^'Hu-)>j] 

of the integrand shows that if u is the point of jump of the process 'D{t), 
and Q(i)(u-) =n<C, then Q^^\u) = 0. 

Similarly to (A.l), equation (II. 2p can be rewritten as follows: 

(A3) Q(2)(t) = A(2)(t)-^ ri{Qa)(^) = o} 

xi|q(2)(^_) >j_ Q(^){u-)\dD{u). 

The explanation of the equivalence of (A.l) and (A. 3) is similar to the above 
case, but slightly more complicated in details. Specifically, the presence of 
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the term I {Q^^'{u) = O} in the integrand is obvious, and the vahdation of 
the term 

l|Q(2)(^_)>j-_g(i)(^_)| 

is explained similarly to that of (A. 2). 

Let us now find the representation for Q2{t) = Q^^'{t) + Q'^-*(t). Keeping 
in mind that Qi{t) = Q^'^^t) and Ai{t) = A^^\t) from (A.l) and (A.3) we 
obtain: 

c 

(u) 



Q2{t) = A2{t) - V / I{Qi(n-) > j}dD{i 

-iZ [ HQiiu-) < j}I{Q2{u-) > j}dD{u) 
= A2it) - V / I{Q2(n-) > j}dD{u). 



The term 



V / I{Qi{u-) > j}dD{u) 
T^ Jo 



characterizes departure lengths from the highest priority buffer, while the 
term 

V / l{Qi{u-)<j]l{Q2{u-)>3]dD{u) 

characterizes that from the second buffer of the lower priority. 
Thus for k = 1,2 we have already shown 

Qk{t) = Ak{t) - V / l{Qk{u-) > j}dD{u). 

Let us prove (j2.ip by using induction. 

For this purpose let us write first a representation for Q*- •*(i), k = 
1,2, . . . ,i — 1. Similarly to (A.3) we have: 

(A4) Q(*^+i)(t) = ^(^'+i)(t)-V f\{Q,{u) = 0} 

xl {Q('^+i)(tx-) > j - Qk{u-)} dD{u). 

Equation (A. 4) is a straightforward extension of (A.3). Therefore, assuming 
that (j2.ip is valid for some k and adding Qk{t) and Q''''^^'{t), and similarly 
to the above for the k + 1st cumulative buffer content we obtain: 

c ,-t 
Qk+i{t) = Ak+iit)-y2 I{Qk+iiu-)>j}dDiu). 



26 VYACHESLAV M. ABRAMOV 

Representation p.ip is proved. 

APPENDIX B: DERIVING (USD 



Equation ()4.4p is a basic equation for our analysis. With cx) • = for 
the left-hand side of (14.41) we have: 



A(*=)(i) 



oo ^ i 

xP|A(^')(t) = /| 

oo -. I 

= A('=)^hmEX;i^l{QW(tf)) >m, Q^'^^ (t^^ 

1=0 i=l 

A('=)(t) 

fc = 1,2, ...,^. 
The right-hand side of (j4.4p can be rewritten as follows: 

C D(t) 
lim -E^^lJQ('=)(xj-) =m-l + H 

1 ^ ft 
lim -EV / i|q('^)(u-) = m-l + /|dD(u), 
^^°o t ^ Jo •- J 



*^^^ - ' < m : 



L, .^, . . . , -I 



t- 
A: = 1,2, 
From these last two equations we obtain: 



S-l4TSi{c'''('f)a-Q'"('f-)<™} 






1=1 ■^^ 
fc = 1, 2, . . . , t, 

and taking into account (II. 6p we finally obtain (j4.5p . 
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